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ABSTRACT 

Aims. We determine Period-Luminosity relations for Milky Way Cepheids in the optical and near-IR bands. These relations can be 
used directly as reference for extra-galactic distance determination to Cepheid populations with solar metallicity, and they form the 
basis for a direct comparison with relations obtained in exactly the same manner for stars in the Magellanic Clouds, presented in an 
accompanying paper In that paper we show that the metallicity effect is very small and consistent with a null effect, particularly in 
the near-IR bands, and we combine here all 111 Cepheids from the Milky Way, the LMC and SMC to form a best relation. 
Methods. We employ the near-IR surface brightness (IRSB) method to determine direct distances to the individual Cepheids after 
we have recalibrated the projection factor using the recent parallax measurements to ten Galactic Cepheids and the constraint that 
Cepheid distances to the LMC should be independent of pulsation period. 

Results. We confirm our earlier finding that the projection factor for converting radial velocity to pulsational velocity depends quite 
steeply on pulsation period, p = 1.550-0.1 86 log(P) in disagrement with recent theoretical predictions. We find PL relations based on 
70 Milky Way fundamental mode Cepheids of Mk = -3.33(±0.09)(log(/') - 1.0) - 5.66(±0.03), Wvi = -3.26(+0.1 l)(log(/') - 1.0) - 
5.96(+0.04). Combining the 70 Cepheids presented here with the results for 41 Magellanic Cloud Cepheids which are presented in an 
accompanying paper, we find = -3.30(±0.06)(log(/') - 1.0) - 5.65(±0.02), Wy, = -3.32(±0.08)(log(P) - 1.0) - 5.92(±0.03). 
Conclusions. We delineate the Cepheid PL relation using 111 Cepheids with direct distances from the IRSB analysis. The relations 
are by construction in agreement with the recent HST parallax distances to Cepheids and slopes are in excellent agreement with the 
slopes of apparent magnitudes versus period observed in the LMC. 

Key words. Stars: variables: Cepheids, Stars: fundamental parameters. Stars: distances, distance scale 

1. Introduction In the present paper we use the new and largely expanded 

] . . data set from Paper II for now 36 LMC Cepheids together with 

, In this series of papers we calibrate the Cepheid Period- jjg^ (jjj.g^.j geometric parallax measurements from Benedict 

> Lummosity (PL-) relation usmg the mfrared surface brightness ^007) to empirically determine the appropriate p-factor 

V ^^^^^^ method. In Paper II we address the effect of metallicity relation to be used in the analyses. 

rS on both the slope and the zero point of the relations in both the „, ^ , 

5-H 1 J Tr. u J J- 11 / • » » We present new accurate radial velocity data for 14 galac- 

fi-i optical and near-IR bands nnding very small (consistent with . „ , . , ,. , , •; ^ „„ 

■ ■ ■ ' zero) effects in the near-IR and small, but possibly significant Cepheids expanding the sample to a total of 77 Cepheids, 70 

effects in the optical which can be used to delineate the Milky Way PL relations 

Gieren et al. (2005) made a first determination of the LMC 
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We have reanalyzed the complete sample using exactly the same 
code and calibrations as for the LMC sample and adopting ex- 
actly the same reddening law to allow a direct comparison. 



PL relations based on thirteen stars with IRSB based distances. 
They found that the distances to the individual Cepheids were 

dependent on the pulsation period which of course is unphysi- Based on the (near-) universality of the PL relations we com- 

cal. They found that the problem could be resolved by changing bine the Milky Way, LMC and SMC samples to determine PL re- 

the adopted projection (/>) factor relation, which converts ob- lotions based on 1 1 1 Cepheids which at the same time constrain 

served radial velocities into pulsation velocities that are needed the slopes very well, and whic h are t ied directly to the parallax 

for Baade-Wesselink type analysis ^^'"o point from Benedict at al. (120071 1. These relations thus form 

a very solid basis for the Cepheid distance scale. 



* Table 3 is only available in electronic form at the CDS The paper is structured as follows: In Sec|2] we present the 

via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via data from the literature as well as new radial velocity data for 
|http: //cdsweb.u -strasbg.fr /cgi-bin/qc at?J/A+A/| fourteen Milky Way Cepheids. In Sec|3] we present the IRSB 



2 



J. Storm et al.: Calibrating the Cepheid Period-Luminosity relation from the infrared surface brightness technique 



method and use the Benedict et al. (120021 12007I I parallaxes as 
well as the results from Paper II on 36 LMC Cepheids with IRSB 
distances to constrain the /7-factor relation. We then use the new 
p-factor relation to determine distances and luminosities for 77 
Cepheids and derive new PL relations for the 70 fundamental 
mode pulsators with good data sets. We proceed to combine the 
data with the Magellanic Cloud sample to give our best global 
PL relations which can be used for distance determination to 
other galaxies. In Sec|4]we compare the results with other recent 
investigations and in Sec|5]we summarize our conclusions. 

2. The Data 

We have searched the literature for optical (V-band) and near- 
infrared (/T-band) light curves as well as radial velocity curves. 
The starting point for the search was the catalogues of complete 
phase coverage /T-band light curves for Milky Way Cepheids 
pubhshed by Welch et al. (fT984l l. Laney & Stobie ( fT992] l and 
Barnes et al. (I1997I I. Since the publication of those papers a rich 
literature of high quality optical and radial velocity data has ma- 
terialized and for the majority of the stars it is now possible to 
apply the near-infrared surface brightness method to determine 
distances and luminosities. 

In addition to the literature data we have obtained new, accu- 
rate, radial velocity curves for 14 of these Cepheids to improve 
the phase coverage and/or data quality for these stars. 

We have selected the data sets according to quality and com- 
pleteness, but also to ensure, as far as possible, that the data have 
been obtained close in time to the near-infrared data to reduce 
possible errors due to period variations of the stars. Known dou- 
ble mode pulsators have been disregarded a priori as the appli- 
cation of the IRSB method to such stars could only be attempted 
if all the data were truly simultaneous. Thus the list of references 
reflects this pre-selection and does not include the data sets that 
were not used for the analysis. A number of first overtone pul- 
sators has also been included, but they are of course not used 
for the delineation of the fundamental mode pulsator Period- 
Luminosity relations. 

In TablUwe present the list of stars and the references to 
the data sets which we have used in the present analysis. A 
more complete list of of data sets can be found in Groenewegen 
(I2OO8I 1. 

The BVIc photometry reported here is all on the Johnson- 
Kron-Cousins system and the different data sets have been trans- 
formed to this system when necessary. Similarly all the near-IR 
data have been transformed to the SAAO system as necessary 
using the transformations from Carter (1 19901 1. 

The radial velocity data have all been obtained using high 
resolution {R > 20000) spectrographs. Most of the radial ve- 
locities have been derived using cross-correlation techniques or 
techniques which are equivalent. In this way the radial velocities 
are assumed to be on a common system and we have not seen in- 
dications of significant systematic diff'erences between datasets 
for any of our stars for which we have had data from different 
techniques. This is an important point for the application of the 
IRSB method as the conversion from radial velocities to pulsa- 
tional velocities, the so called /^-factor, depends to some extent 
on the procedure which was used for deriving the radial velocity 
from the input data. 

2.1. New radial velocity measurements 

For fourteen of the stars the radial velocity curves were either 
missing or they were of limited quality. We have obtained 381 



Table 2. The Cepheids for which we have obtained new radial 
velocity curves. For each star the number of new data points is 
listed as well. 



Star 


A'obs 


Star 


A'obs 


TT Aql 


25 


LS Pup 


26 


V496 Aql 


19 


VZ Pup 


29 


4- Gem 


65 


XSgr 


19 


VZ Cyg 


20 


YSgr 


15 


XPup 


42 


BB Sgr 


18 


AQPup 


38 


XXSgr 


20 


BNPup 


31 


YZSgr 


14 



Table 3. The heliocentric radial velocities (RV) measured with 
the STELLA echelle Spectrograph (SES) as returned by the data 
reduction pipeline and offset to the CORAVEL velocity zero 
point. The complete table is available in electronic form from 
the CDS. 



Star 


HJD 


Phase 


RV 


cr(RV) 




(Days) 




(kms"' ) 


(kms"' ) 


TTAql 


2454175.76609 


0.53 


14.79 


0.04 


TTAql 


2454213.73669 


0.29 


-8.15 


0.06 


TTAql 


2454218.64326 


0.65 


24.95 


0.07 


TTAql 


2454222.59739 


0.94 


3.25 


0.02 


TTAql 


2454223.59473 


0.01 


-6.92 


0.02 



new radial velocity measurements for these stars (see Tab|2]) 
using the STELLA Echelle Spectrograph (SES) mounted on 
the fully robotic 1.2m STELLA-I telescope (Strassmeier et al. 
2004 2010 and Weber et al. 2008) at the Izana Observatory on 
Tenerife, Spain. SES is a fiber-fed echelle spectrograph with a 
2kx2k CCD detector covering the wavelength range from 388 
to 882 nm with small inter-order gaps starting at 732 nm and 
increasing towards the red. The resolving power is R = 55000 
giving a spectral resolution of 0.12 A at 650 nm. 

The spectra were obtained in fully robotic mode in the period 
from March 2007 until July 2010 and reduced using the auto- 
matic data reduction pipeline (Ritter & Wasshiittl, 120041 Weber 
et al. l2008l l developed for the instrument. 

The radial velocities returned by the pipeline were corrected 
for instrumental velocity offsets and placed on the CORAVEL 
system by applying the offset of -i-0. 503km s"' determined by 
Strassmeier et al. (in prep). The radial velocities are tabulated in 
Tabl3]and plotted in FiglT] 

2.2. Pulsation velocities 

The p-factor (see Sec l3.4b . which is used to convert the observed 
radial velocities into pulsational velocities, depends to some ex- 
tent on the spectrograph and the procedure used for extracting 
the velocities as different spectral features might carry different 
weight in deriving the pulsational velocity. To investigate this ef- 
fect for the STELLA velocities we have observed the star TT Aql 
for which an excellent CORAVEL based velocity curve is avail- 
able (Imbert 1999). 

We follow the procedure described by Storm et al. (2004| 
to determine the ratio between the p-factors for STELLA and 
CORAVEL based velocities, //stella = Pstella/pcor, for the 
relevant phase interval, < < 0.8. In Fig|2]we plot the differ- 
ence in pulsational velocity, 

AVpuis = /?cor(RVcor - Vy) - 77/'cor(RVstella -Vy) (1) 
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ID 


Optical 


A^-band 


rad. vel. 


ID 


Optical 


A^-band 


rad. vel. 


77 Aql 


4,17 


3,4 


9,12,29,35 


UNor 


7,16,20 


2 


8,21 


U Aql 


17 


3 


32 


QZ Nor 


6,47 


2 


28,41 


FF Aql 


13,17 


3 


53 


TW Nor 


16,20 


2,3,58 


23,41 


FM Aql 


4,17 


3,4 


35 


V340 Nor 


42 


2 


23 


FN Aql 


17,20 


3 


34 


Y Oph 


7,17,20 


2 


27,55 


SZ Aql 


17,20 


2 


35,40 


BF Oph 


10,17 


2 


11,14,15 


TT Aql 


4,7,13,17,20 


3,4 


1 


XPup 


17,20 


2 


1 


V496 Aql 


10,17,47 


3 


1,48 


AQPup 


7,17 


2,22,24 


1 


RT Aur 


13,17 


4 


12,34 


BN Pup 


7,16,20 


2,22 


1 


f Car 


20,38 


2 


40 


LS Pup 


7 


2,22 


1 


UCar 


7,20 


2,3 


8,40 


RS Pup 


17,20,47 


2 


9 


V Car 


20 


2 


15 


VZPup 


7,20 


2,22 


1 


CT Car 


51 


2,22 


21,41,47 


KQ Sco 


7 


2,3 


8,54 


VY Car 


7 


2 


8,15 


RY Sco 


7,16,17,20 


2 


8,14,15 


WZ Car 


7,20 


2,22 


8 


EV Set 


17,20 


2 


9,19,23 


SU Cas 


4,17 


4 


9,23 


RU Set 


17,20 


2,3 


19,41 


KNCen 


7 


2 


8,21 


SS Set 


17,20 


3 


54 


VCen 


10,20 


2 


11,15 


SSge 


4,13,17 


3,4 


34 


VW Cen 


7,20 


2 


8 


GYSge 


43 


2 


19,34 


XXCen 


39 


2 


8 


USgr 


10,17 


2 


9,23 


(5Cep 


4,17 


4 


9,12,23,29,35 


W Sgr 


17,20 


3,59,60 


23 


SU Cru 


20 


2 


8,21 


XSgr 


17,20 


3,57 


1 


XCyg 


13,17 


3,4 


9,29 


YSgr 


17,20 


3 


1,56 


CD Cyg 


17 


3 


26 


BB Sgr 


10,17 


2 


1 


DTCyg 


17 


3 


23 


WZSgr 


17 


2,3 


34 


SU Cyg 


13,17 


3 


31 


XX Sgr 


17 


3 


1 


VZ Cyg 


4,17 


3,4 


1 


YZSgr 


17,47 


3 


1,48 


ySDor 


20,33 


2 


37 


V350 Sgr 


10,17 


3 


11,27,34 


f Gem 


13 


57 


1 


SZTau 


4,17 


2,4 


23 


XLac 


4,17 


4 


23 


T Vel 


47,20 


2 


11,15 


YLac 


4,17 


4 


26 


CS Vel 


44,45 


2,22 


23,41 


J_,ac 


A 1 7 
4, 1 / 


A 
4 


JO 


K I vei 


/,lo 




o,Zi 


BGLac 


4,17 


4 


26 


RZ Vel 


7,16,20 


2 


8,56 


GHLup 


7,20 


2 


8 


SW Vel 


7,16,20 


2 


21,27 


TMon 


47,7 


2 


9,23 


S Vul 


43 


2,3 


34 


CV Mon 


17,20 


2 


9,41 


T Vul 


13,17 


3,4 


23 


SMus 


46 


2,3 


30 


U Vul 


4,13,17 


4 


23,25 


UU Mus 


16,39 


2 


8 


SV Vul 


13,17 


2,3 


23,26 


SNor 


44,46 


2 


23 











(1) This paper; (2) Laney & Stobie. , 1992, (3) Welch, et al. ri984l (4) Bames et al. ri997l (6) Berdnikov & Tumer[T995] (7) Coulson, & Caldwell[T985] (8) Coulson, Caldwell, & Gieien 
[19851 (9) Storm et al. l2004l (10) Gieren [T98Tbl (11) Gieren [T98Tal (12) Ki5s ll998bl (13) Kiss 1998a (14) Barnes, Moffett, and Slovak fMSl (15) Lloyd Evans. fT980l (16) Madore fT975l 
(17) Moftett & BarnesfTgOT (19) Metzger et al.[T99T] (20) Pel [T976l (21) Pont, Mayor, & Burki l 19941 (22) Schechter et al. [T992l (23) Bersier et al. [T994l (24) Welch [T985l (25) Imbert 
fT996l (26) Imbert 1999 (27) Petterson et al. l2005l (28) Kienzle et al.. fT999l (29) Butler & Bell fT997l (30) Evans [T990l (31) lmbert fT984l (32) Welch et al. flWTl (33) Shobbrook|T992] 
(34) Gorynya et al. il998. (35) Barnes et al. l2005l (36) Sugars & EvansTW (37) Wallerstein et al^WT (38) Bersier l2002l (39) Coulson, Caldwell, & Gieren [T985l (40) Bersier l2002l 

(41) Metzger et al. [T992l (42) Bersier et al. [T994l (43) Berdnikov 1986 1987, 1992, 1992 1992, 1992, 1992] (44) Berdnikov & Tumer [T998l (45) Berdnikov & Turner l2000l (46) 
Walraven, Tinbergen, & Wakaven fT964l (47) Berdnikov & Caldwell 2001 (48) Caldwell et al. 2001 (51) Pojmanski et al.; (53) Evans fT990l (54) Groenewegen'2008' (55) Nardetto et al. 
l2006l (56) Nardetto et al. l2009l (57) Feast et al. [2008, (58) McGonegal et al. ,1983, (59) Kimeswenger et al. l2004l (60) Wisniewski & Johnson ,1968. 



as a function of phase between the linearly interpolated observed 
radial velocities for the two spectrographs for three different val- 
ues of 7]. It appears that in the phase interval from 0.15 to 0.75 
where the velocity difference curve is smooth, the formally best 
value is 77 = 1.03. It is however also clear that the uncertainty 
is rather large and given that the datasets have been obtained at 
epochs differing by about ten years, we choose not to apply any 
additional corrections to the STELLA velocities but assume that 
the STELLA and CORAVEL p-factors agree to within 3%. For 
the present sample of 14 stars with STELLA velocities a 3% ef- 
fect on the p-factor translates into a 1 % effect on the slopes of 
the final PL relations which is much smaller than the statistical 
errors. 

For some stars we have used the radial velocity data from 
Nardetto et al. (2009) using the HARPS data pipeline. This 
data set also contains data for a star, ^ Gem, for which a 



good CORAVEL data set is available from Bersier ('1994JI. 
Unfortunately the period of the star is not perfectly constant 
over time so it has been necessary to shift the two radial ve- 
locity curves with respect to each other to obtain a good match. 
Proceeding then as for the STELLA dataset we find that the two 
data sets are in good agreement and that t/harps - 1 -00 + 0.03 as 
can be seen in Figl3] 



For another star, f Car, there is a similar possibility of a direct 
comparison between CORAVEL and HARPS data. The radial 
velocity curves from different data sets exhibit some variations 
though, and the conclusions are less straight forward than for 
f Gem, but they do agree with a value of t] which is indistin- 
guishable from unity, which we will adopt in the following. 
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Fig. 1. The new radial velocity curves for fourteen Milky Way Cepheids from the STELLA echelle spectrograph as tabulated in 
Tabs 




Phase Phase 

Fig. 2. The difference in pulsational velocity as a function Fig. 3. The difference in pulsational velocity as a function of 
of phase between CORAVEL and STELLA measurements for phase between CORAVEL and HARPS measurements for three 
three different choices of 77 where // - pstella/ /'cor- different choices of 77 where 77 - Pharps/pcor- 

3. The Analysis 

3.1. The IRSB method 

The infrared surface-brightness (IRSB) method is a variant of the 
Baade-Wesselink method originally developed by Barnes and 
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Evans (I1976I I in the optical wavelengths. It is based on a func- 
tional relation between a colour index and the surface bright- 
ness parameter in the V-band, Fy- It was originally calibrated by 
Welch ([1994) and a few years later Fouque and Gieren ( I1997I I. 
They determined a very tight linear relationship between the 
{V-K) colour index and Fy based on interferometric angular di- 
ameters of giant stars found in the literature, thus extending the 
method to the near-infrared. The scatter in this relation was sig- 
nificantly smaller than was the case for the optical colour indices 
used previously. A detailed description of the implementation of 
the method which we use in the present paper can be found in 
Storm et al. (2004|. 

Recently direct interferometric angular diameter measure- 
ments of Cepheids have become available (Nordgren et al. 120021 
Kervella et al. l2004al Merand et al. [2005i ) allowing a direct com- 
parison between the surface-brightness relation for static stars 
with actual pulsating stars. On this basis Kervella et al. (2004b) 
find excellent agreement between static and pulsating stars as 
well as with the relation determined by Fouque and Gieren 
d 19971 ) for static stars. They find a best fit relation of 



-0.1336(y-/i:)o + 3.9530 



(2) 



with the coefficients determined to better than 2%. We adopt 
their relation for the following analysis. 

The surface brightness measure Fy is directly related to the 
stellar angular diameter, 6 through the relation 



Fv{^) = 4.2207 - 0. 1 yo(0) - 0.5 log 9{(l>) 



(3) 



where Vq is the de-reddened visual magnitude, and is the 
phase. 

At the same time geometry gives us the stellar angular diam- 
eter from the stellar radius, Risf)), through the relation 



e{(p) = 2R{(l))ld = 2{Rq + AR{<p))/d 



(4) 



where tp is the pulsation phase, d is the distance and R is the 
radius. 

Integrating the radial velocity curve then gives the radius 
variation, ARicf)) between a reference radius, Rq and the given 
phase, as 



AR(cf>) 



(5) 



where p is the so called projection factor converting radial veloc- 
ity into pulsational velocity, V, {<p) is the observed radial velocity 
and Vy is the systemic velocity. 

We can now solve Eq|4]for the two parameters, mean radius, 
and distance, d by linear regression to the observed values of 
0(0) from Eql3] and AR((p) from Eq|5] As discussed in Storm et 
al. (2004) we use the OLS bi-sector fit from Isobe et al. ( 199()| 
for the regression fit. We fit only the phase interval (p e [0.0, 0.8] 
where the shapes of the two curves usually agrees very well, 
and we avoid the phase region (/> G [0.8, 1.0] where the agree- 
ment often is quite poor, most likely due to shocks in the stellar 
atmosphere. We also allow for a small phase shift between the 
photometric and radial velocity data to optimize the quality of 
the fit. The effect on the final PL relation of these phase shifts is 
mainly to decrease the scatter in the relation. 

An example of the fit for the star BB Sgr is shown in FiglH In 
the upper panel the data used for the actual OLS bi-sector fit (see 
Storm et al. l2004l for more details on this) can be seen and in the 
lower panel the corresponding photometric angular diameters 




0.4 0.6 
Phase 



Fig. 4. The IRSB fit to the data for the star BB Sgr. The deviation 
between photometric and spectroscopic angular diameters in the 
phase interval </> e [0.8, LO] is evident and is observed for many 
of the Cepheids in our sample. This phase interval is therefore 
disregarded in the fit for all the stars. 



have been plotted as filled squares for the points used in the fit 
and red crosses for the points in the phase interval 6 [0.8, 1.0] 
which have been disregarded in the fits. The curve in the lower 
panel delineates the corresponding spectroscopic angular diam- 
eter. 



3.2. Absorption 

In order to derive dereddened magnitudes for our Cepheids we 
use the colour excess values as given in Fouque et al. (2007}. 
These values are on the system defined by Laney & Caldwell 
( 120071) and as discussed by Fouque et al. ( 120071) . these values 
are in agreement with the system established by Tammann et al. 
(120031) . who recalibrated the original measurements compiled by 
Fernie ^995 ). 

For the reddening law we similarly adopt the choice made 
by Fouque et al. ( 120071 ). namely the law determined by Cardelli 
et al. ( 119891) with a total-to-selective absorption in the V band of 
Ry - 3.23 as determined by Sandage et al. (120041) . For the other 
bands we use A/ = 0.6QSxAv,Ak = 0.119xAv,Ay = 0.292xAv, 
and Ah = 0.181 xAy. 



3.3. Fourier coefficients and identification of fundamental 
mode pulsators 

Before we attempt to establish the PL relations we have to iden- 
tify the fundamental mode pulsators in the sample. This is diffi- 
cult without referring to the period-luminosity diagram. We use 
the Fourier parameters for the radial velocity data and the di- 
agrams from Kienzle et al. (T999\ to reject overtone pulsators. 
To be conservative we use the A] parameter (the radial velocity 
amplitude) to reject .s-type Cepheids from our final sample as 
well. In this way we might remove some bona fide fundamental 
pulsators as well, but we ensure that we have a uniform sample. 
The Fourier parameters for all the Cepheids based on the data 
sets given in TablU are tabulated in Tab|4] The eiTor estimates 
on the parameters have been calculated using the approximative 
formula given by Petersen (I1986I ). 
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Fig. 5. The Fourier parameters based on the radial velocity data 
for the Cepheids with periods less than 9 days. The dashed line 
indicates the adopted division between normal Cepheids and s- 
Cepheids based on Fig. 2 in Kienzle et al. ( 1999). Filled circles 
indicate fundamental pulsators and the triangles, first overtone 
pulsators. 

Table 5. The distances for Cepheids with HST parallax mea- 
surements from Benedict et al. ( I2007I I. The HST distances are 
given in column 3, the IRSB distances using our preferred p- 
factor relation from Eq|7] are given in column 4 and the differ- 
ence in column 5. The uncertainty on the difference is given in 
the last column. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


ID 


log(P) 


rf(HST) 


^/(IRSB) 


Ad 


cr(Arf) 






(pc) 


(pc) 


(pc) 


(pc) 


FF Aql 


0.650390 


355.9 


369.8 


14.0 


42.3 


RT Aur 


0.571489 


416.7 


389.0 


-27.6 


33.0 


^ Cai- 


1.550820 


497.5 


517.6 


20.1 


47.9 


(5Cep 


0.729678 


273.2 


266.7 


-6.5 


20.3 


/3Dor 


0.993131 


318.5 


326.6 


8.1 


22.0 


^Gem 


1.006500 


359.7 


385.5 


25.8 


38.3 


W Sgr' 


0.880529 


438.6 


216.8 


-221.8 


71.1 


XSgr 


0.845893 


333.3 


322.1 


-11.2 


28.1 


YSgr 


0.761428 


469.5 


436.5 


-33.0 


75.6 


T Vul 


0.646934 


526.3 


542.0 


15.7 


59.4 



t Not considered in the fits. This star is a known spectroscopic binary. 

In Fig|5]the Fourier parameters for the short period Cepheids 
have been plotted and the i-Cepheids identified as the stars with 
Al below the dashed line. These stars are DT Cyg, EV Set, 
FF Aql, SU Cas, SZ Tau, and QZ Nor. SU Cas is the shortest- 
period Cepheid in our sample and was probably the first Galactic 
Cepheid for which pulsation in the first overtone mode was 
firmly established by Gieren j 19761 [1982 ). Two stars lie below 
the sequence of fundamental mode pulsators in the Ai vs. P plot 
but above the dashed line and in the other Fourier parameters 
they appear unremarkable, so they are not obvious i-Cepheids 
or overtone pulsators. These stars are X Lac, and V496 Aql and 
we will keep them in the sample of fundamental mode stars. 

3.4. Constraining tlie projection factor 

For any Baade-Wesselink type method it is necessary to deter- 
mine the pulsational velocity of the surface of the star as this is 



T 




HST parallax zero-point 



1.2 1.4 1.6 1.8 

Fig. 6. Constraints on the /j-factor relation from the HST par- 
allax of Milky Way Cepheids (full line), and from requiring the 
LMC distance to be independent of the pulsation period of the 
stars (dashed line), see Sec l3.4l for details. The one sigma con- 
tours are also shown with thin lines. The filled square with eiTor 
bars shows the theoretical relation from Nardetto et al. (I2009I I. 
the open triangle the Hindsley and Bell (1 19861 1 relation and the 
filled circle the best fit. The open square shows the theoretical 
constraint shifted in /3 to comply with the HST parallax values. 

the velocity curve which is matched against the angular diameter 
curve from the photometry. The same is the case when the an- 
gular diameter curve is directly determined from interferometry. 
The conversion from the observed radial velocity to pulsational 
velocity is commonly parametrized with the projection factor, p. 

The /7-factor is largely a geometrical coiTection taking into 
account the fact that the radial velocity that we measure is based 
on light coming from all points of the hemisphere of the star 
which is visible to the observer and not just from the surface 
element which move along the Une of sight to the star. In fact 
as discussed by Sabbey et al. d 19951 1 the p-factor depends on the 
temperature of the star, which changes with pulsation period, as 
it depends on the limb darkening of the star. 

The /9-factor used for converting the observed radial veloc- 
ities into pulsational velocities has direct consequences for the 
derived distances, as it scales directly with the stellar radius vari- 
ation, and is probably the largest source of systematic error for 
the method. In the past we (e.g. Gieren et al. 119931 Storm et al. 
i200 4, Barnes et al. .2005 ) have used a relation with a weak period 
dependence, p = 1 .39 - 0.03 log P, based on theoretical work by 
Hindslev and Bell ( fT986] l. 

Later we (Gieren et al. 120051 ) found that the use of this p- 
factor relation with the IRSB method for LMC Cepheids leads 
to an unphysical dependence of the distance modulus on the 
pulsation period. We found that a /j-factor relation of p - 
1.58 - 0.15 logP removed this period effect, but the conclusion 
was not very firm due to the limited sample of only 13 LMC 
Cepheids. 

On the theoretical side, Nardetto et al. (120071 ) have care- 
fully analyzed individual line profiles from pulsating atmosphere 
models and they also found a relation which was steeper than the 
Hindsley and Bell (1986) slope. Nardetto et al. (2009) general- 
ized this work to match the cross-correlation technique used in 
most observational work on radial velocities and found a relation 

p = 1.31(+0.06)-0.08(+0.05)logP (6) 
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With the sample of 36 LMC Cepheids, presented in Storm 
et al. (1201 II hereinafter Paper II), covering a wide range of pe- 
riods we are in a much better position to constrain the /^-factor 
relation empirically. Furthermore the recent direct parallax mea- 
surements for ten MiUcy Way Cepheids by Benedict et al. (120021 
I2007I I using the Hubble Space Telescope Fine Guidance Sensor 
provide a fundamental set of reference data which we can use to 
calibrate the /^-factor relation. 

For three of the ten stars with HST parallaxes (Y Sgr, X Sgr, 
and ^ Gem) we present new radial velocities here, significantly 
improving the available data quality. For one star, W Sgr, which 
is a known binary (Szabados 2003) the IRSB fit is clearly very 
poor and we disregard this star in the further analysis. X Sgr is 
known to be affected by a double Shockwave in its atmosphere 
(Mathias et al. 2006) but the IRSB fit looks fine and the agree- 
ment with the HST parallax is excellent so we keep it in the 
sample. This leaves us with nine stars in common. 

As a first step we use the theoretical relation from Nardetto 
et al. (I2009I I in Eq|6] and apply the IRSB method to these nine 
stars. We find a disappointing difference of -0.30 + 0.05 mag in 
the distance moduli, the IRSB distances being shorter. Applying 
the Nardetto et al. relation to the LMC Cepheids in Paper II we 
similarly find an unlikely result, namely an LMC distance mod- 
ulus of 18.26 + 0.04, much shorter than the canonical value of 
18.50. So if we proceed using first principles, we have a seri- 
ous conflict with the Benedict et al. (2007) result as well as with 
most recent works on the LMC distance that have confined the 
true distance modulus to a value between 18.4 and 18.6 (e.g. 
Pietrzyiiski et al. 2009; Szewczyk et al. l2008] l. 

To reconcile these results we have to conclude either that the 
theoretical /^-factor relation is incorrect due to the lack of some 
physics, or that there is a period dependence in the IRSB method 
which is not properly accounted for in the current theoretical p- 
factor relation. We attempt to constrain this effect empirically 
and parametrize it as a part of the p-factor relation which we 
retain as having the simple linear form p - a^jX log(P) + j3p. 

We need to determine two parameters, namely the slope (ap) 
and the zero point (J3p) of the p-factor relation. We have two 
independent constraints, namely that there should be no system- 
atic dependence of the LMC Cepheid distances with pulsation 
period, and we should reproduce, on average, the Benedict et al. 
( I2007I I distances. 

It turns out that these two constraints are largely orthogo- 
nal in the (J3p, Op) plane as can be seen in Fig|6] To determine 
these parameters we simply apply the IRSB method to each 
sample of stars (Milky Way and LMC) for an array of slopes 
and zero points for the /j-factor relation and see where the con- 
straints are fulfilled. We have varied the zero-point in the range 
j3p € [1.2, 1.75] in steps of 0.05 and the slope of the relation in 
the range Up e [-0.31, 0.01] in steps of 0.03. 

In Tab|5]we have listed the stars with parallax distances from 
Benedict et al. (2002 2007). For each pair {ap,/3p) we have deter- 
mined the IRSB distance, (i(IRSB), to these stars and computed 
the difference A(d) = t/(HST) - <i(IRSB). To weight the points 
independently of distance, we have normalized the values by di- 
viding by the average distance, davg = (li(HST) + c/(IRSB))/2., 
before computing the offset A{d)ldavg- We then took the mean 
value of these offsets and determined the values in the fip - ap 
plane where this mean offset is zero. This is a straight line which 
is shown in Fig|6] as a full line with the two thin parallel lines 
showing the estimated 1-cr interval. 

We then turn to the LMC data set and proceed as for the 
Milky Way sample and carry out the IRSB analysis for the same 
set of {ap,/3p) values. We then look for the points where the slope 



of the LMC distance modulus as a function of log(f ) is zero. For 
each individual Cepheid distance we apply the distance modulus 
correction A(m - M) from van der Marel and Cioni (2001) to 
correct for the inclination of the LMC disk before determining 
the slope. In the /3p, ap plane the resulting constraint is shown as 
a dashed line and the two thin dashed lines indicate the estimated 
1-cr interval. 

The best estimate is: 

p = 1.550(+0.04) - 0.186(+0.06)logP (7) 

This is shown as a filled circle in Fig|6] The /9-factor law is even 
a bit steeper with period than the relation which we found earlier 
(Gieren et al. 2005) and it differs even more from the recent the- 
oretical relation from Nardetto et al. (2009) which is shown as a 
filled square in Fig|6] For reference the relation from Hindsley 
and Bell (1986) has been plotted as an open triangle in the fig- 
ure. It agrees within one cr with the HST parallax constraint but 
not with the constraints from the period dependence for the LMC 
sample. Additionally, the open square shows a relation where we 
have adopted the slope ap = -0.08 from Nardetto et al. (2005|l 
but forced the zero point, /3p, to give agreement with the HST 
parallaxes. 

We can look at the problem in a slightly different way and 
determine the /7-factor for each of the Cepheids with measured 
HST parallax distances, and plot them as a function of log(P) by 
forcing the IRSB distance to be equal to the parallax distance. 
We have done this and show the results in Fig|7] A linear fit 
to these values gives p = -0.28(+0.08) log(P) + 1.65(+0.07). 
Within the errors this agrees with our relation (Eq|7]i. We pre- 
fer, however, to use the relation in Eq|7] as it is based on many 
more stars, especially at pulsation periods longer than ten days. 
The linear fit is shown in the figure as well as our adopted re- 
lation (labelled LMC and HST ZP) and the Hindsley and Bell 
(1986) and Nardetto et al. (2009) relations. The horizontal line 
atp- 1 .5 indicates the limit above which the /^-factor would im- 
ply an unphysical limb-brightening instead of the expected limb- 
darkening. We note that for the short period stars the p-factor is 
coming close to this limit. 

In Tab|6] we have summarized the resulting values of true 
LMC modulus, the distance offset to the HST parallax data, the 
slope of the LMC Cepheid moduli as a function of period, for 
the different assumed p-factor relations. From this Table we can 
see that the Hindsley and Bell (1986) relation leads to more than 
2cr deviation for the slope of the LMC distance modulus and 
thus seems to be ruled out. The Nardetto et al. ( I2009I I relation 
disagrees on both the constraints, and seems to be ruled out as 
well. Changing the value of /3p to /3p - 1 .455 brings the distance 
zero point into agreement with the HST parallax value, but still 
the distances to the LMC Cepheids are significantly dependent 
on the pulsation period. We thus adopt the fitted relation from 
Eq|7]in the following. 

In Figl8]we have plotted the distance difference between the 
HST parallax distances and our IRSB based distance when us- 
ing the revised p-factor relation. We note that the scatter is very 
small and that the data are consistent with no period dependence 
of the differences. 

3.5. The Period-Luminosity relations 

Using the p-factor relation derived in the previous section to- 
gether with the reddenings discussed in Sec l3.2l we obtain the 
distances and absolute magnitudes for our Milky Way Cepheids 
as given in Tab|7] In that Table we present the adopted pulsation 
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Fig. 9. The Period-Luminosity relations based on the IRSB luminosities for our sample of Milky Way Cepheids in the K- and 
V-bands as well as in the Wesenheits indices. The filled circles represent the fundamental mode Cepheids, and the crosses indicate 
stars which have been disregarded in the linear regression for reasons mentioned in the text. 



period as well as the distance modulus with the formal uncer- 
tainty from our OLS bi-sector fit (see Storm et al. 120041 for de- 
tails). As discussed by Barnes et al. ( 12005l l these uncertainties 
are underestimated by on average a factor of 3.4 when compared 
to the uncertainties returned by the Bayesian fitting technique 
employed in that paper. In the columns 7 to 12 we give the ab- 
solute magnitudes in the B, V, I, J, H, &K bands and in column 
13 and 14 we give the Wesenheit indices (Madore 1 19821 1 in the 
(V - I) band defined as Wvi = My - 2.54(Mv - M/ ) and in 
the near-IR J and K band, Wjk ^ Mk - 0.688(M., - Mk ). In 
column 15 we give the adopted reddening value and in column 
16 the adopted phase shift between the radial velocity data and 
the photometric data (See Storm et al. 2004 for details). 

Using the data from Tab|7] for the fundamental mode pul- 
sators, and after eliminating the binary Cepheid W Sgr which 
exhibits a very poor IRSB fit, we are ready to determine the PL 
relations in the different bands. In Tab[8]we list the resulting re- 
lations together with the observed dispersions around the fits. 
These are our best estimates of the Milky Way Cepheid PL rela- 
tions. In Fig|9]we plot the PL relations in the K, and V bands as 
well as in the Wesenheit indices Wvi and Wjk ■ We note that the 
dispersion around the fits range between 0.22 and 0.39 mag, the 
B-band relation showing a significantly larger dispersion than 
the others. For the other bands the dispersion is only weakly 



Table 8. The Period-Luminosity relations for the Milky Way 
Cepheids of the form M - ax (log(P) - 1 .0) -i- b. The dispersion 
around the fit is also tabulated as well as the formal uncertainties 
on the coefficients as returned from the linear regression. 



Band 


a 




b 


cr 


Mb 


-2.13 ± 


0.13 


-3.28 + 0.05 


0.39 


Mv 


-2.67 ± 


0.10 


-3.96 + 0.03 


0.26 


Ml 


-2.81 ± 


0.10 


-4.76 + 0.03 


0.23 


Wv, 


-3.26 ± 


0.11 


-5.96 + 0.04 


0.26 


Mj 


-3.18 ± 


0.09 


-5.22 + 0.03 


0.22 


Mh 


-3.30 ±0.08 


-5.59 ±0.03 


0.22 


Mk 


-3.33 ± 


0.09 


-5.66 ± 0.03 


0.22 


Wjk 


-3.44 ± 


0.09 


-5.96 ± 0.03 


0.23 



dependent on the wavelength suggesting that the errors on the 
absolute magnitudes are dominated by the distance errors rather 
than by the intrinsic width of the PL relation and/or errors in the 
absorption coiTections. We note that Persson et al. (2004) found 
a dispersion of only 0.11 mag for the /T-band relation in the LMC 
whereas in Paper II we have obtained a value of 0.22 mag, again 
suggesting that the dispersion in our current work is dominated 
by distance errors rather than intrinsic luminosity variations be- 
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Fig. 7. The /j-factor values derived for each of the HST parallax 
Cepheids with the only constraint being that the IRSB distance 
should agree with the parallax distance. The full line shows the 
linear regression fit to the points. The line labelled "LMC and 
HST ZP" shows our adopted relation based on the constraints 
from the LMC Cepheids and forcing the distance zero point 
to agree with the HST parallax values. The stippled line de- 
lineates the theoretical relation from Nardetto et al. (I2009I I and 
the long dashed line the classical Hindsley and Bell ( 119861) re- 
lation. The horizontal line at p - 1 .5 shows the physical limit 
above which the theoretical /j-factor would indicate an unphysi- 
cal limb-brightening. 
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log(P) 



Fig. 10. The Period-Luminosity relation in the /T-band for the 
complete sample of Milky Way, LMC and SMC Cepheids hav- 
ing IRSB-determined distances in our papers. The overtone pul- 
sators and stars which have been eliminated for other reasons as 
discussed earlier have been eliminated from the plot for clarity. 

Table 9. The Period-Luminosity relations based on the com- 
bined samples of Milky Way, LMC and SMC Cepheids in the 
form M - fl(log(f) - 1.0) H- b. The metallicity effect, y, on 
the zero point determined in Paper II is given in the last col- 
umn. The estimated uncertainty on the y value is estimated to be 
0.10 mag dex ' in Paper II. 



100 - 




-1 
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log(P) 

Fig. 8. The difference between the geometrical parallax distance 
measures for nine Cepheids from Benedict et al. (l2007i i and the 
IRSB distances values, plotted as a function of log P. The re- 
vised p-factor relation determined in this paper has been used to 
calculate the IRSB distances to the stars. 



tween Cepheids of similar periods in the sample due to the finite 
width of the instability strip, or eiTors in the reddenings. 

3.6. The combined sample 

In Paper II we show that the PL relations for the Milky Way and 
Large Magellanic Cloud samples are identical within the uncer- 
tainties, particularly in the near-IR bands. This means that we 



Band 


a 


b 


Std.dev. 


r 




(mag dex"' ) 


(mag) 


(mag) 


(mag dex"' ) 


Mv 


-2.73 ± 0.07 


-3.97 ± 0.03 


0.26 


+0.09 


M, 


-2.91 ± 0.07 


-4.75 ± 0.02 


0.23 


-0.06 


Wv, 


-3.32 ±0.08 


-5.92 ± 0.03 


0.26 


-0.23 


Mj 


-3.19 ±0.06 


-5.20 ± 0.02 


0.22 


-0.10 


Mk 


-3.30 ± 0.06 


-5.65 ± 0.02 


0.22 


-0.11 




-3.38 ±0.06 


-5.96 ± 0.02 


0.23 


-0.10 



can combine the data from the two papers to derive a PL relation 
based on a total of 1 11 Cepheids. In Fig [10] we have plotted the 
K-hand absolute magnitudes for the Milky Way, LMC and SMC 
Cepheids from the two papers together. The agreement is excel- 
lent, and a linear regression to the combined sample leads to a 
best determination for the /T-band PL relation of: 

Mk = -3.30(+0.06)(logf - 1.0) - 5.65(+0.02) (8) 

with a dispersion of 0.22 mag. Due to the limited metallicity de- 
pendence of this relation found in Paper II, this relation can be 
directly used for distance determination to galaxies with metal- 
licities between SMC and solar. 

We have listed the combined relations in the other bands, 
including the Wesenheit indices in Tab|9] in all cases without 
applying any metallicity corrections to the absolute magnitudes. 
The relations have been used in Paper II to determine the PL 
relation zero-point dependence on metallicity, y, and for conve- 
nience we have tabulated those values here as well. 
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Table 4. The Fourier parameters for the stars based on the radial velocity data. 



ID 


Period 

(day) 


Systemic velocity 

(kms"' ) 


Al 

(kms-') 


R21 




021 


<hi 


o-(J?2i) 


0-(R3l) 


0-(<*2l) 


0"((*3l) 


SU Cas 


1.94932 


-7.36 


8.90 


0.220 


0.053 


4.348 


2.576 


0.006 


0.006 


0.157 


0.235 


DTCyg 


2.49921 


-1.70 


6.93 


0.203 


0.045 


4.434 


2.067 


0.012 


0.012 


0.218 


0.326 


EV Set 


3.09099 


17.26 


7.50 


0.185 


0.046 


4.741 


2.748 


0.012 


0.012 


0.221 


0.330 


SZTau 


3.14895 


-0.65 


9.42 


0.259 


0.075 


5.024 


3.784 


0.013 


0.013 


0.233 


0.346 


SS Set 


3.67125 


-7.24 


14.90 


0.298 


0.097 


4.449 


2.616 


0.002 


0.002 


0.096 


0.143 


RTAur 


3.72832 


19.57 


15.36 


0.319 


0.145 


4.618 


3.057 


0.004 


0.004 


0.131 


0.194 


QZNor 


3.78655 


-39.50 


7.33 


0.176 


0.043 


5.186 


4.732 


0.012 


0.012 


0.223 


1.855 


SUCyg 


3.84549 


0.08 


16.38 


0.334 


0.147 


4.545 


2.836 


0.010 


0.010 


0.209 


0.303 


BF Oph 


4.06751 


-28.90 


14.81 


0.222 


0.137 


4.420 


3.110 


0.035 


0.035 


0.423 


0.580 


YLac 


4.32378 


-22.20 


16.16 


0.340 


0.192 


4.570 


2.884 


0.010 


0.010 


0.208 


0.303 


T Vul 


4.43541 


-2.90 


14.32 


0.358 


0.162 


4.517 


2.829 


0.004 


0.004 


0.123 


0.185 


FF Aql 


4.47085 


-0.09 


7.83 


0.129 


0.057 


5.108 


3.615 


0.039 


0.039 


0.400 


0.593 


T Vel 


4.63982 


5.75 


14.72 


0.314 


0.103 


4.665 


3.037 


0.038 


0.038 


0.446 


0.615 


VZCyg 


4.86445 


-18.30 


15.39 


0.392 


0.155 


4.716 


3.302 


0.002 


0.002 


0.082 


0.123 


V350 Sgr 


5.15424 


7.73 


15.65 


0.380 


0.193 


4.701 


3.320 


0.004 


0.004 


0.129 


0.191 


BGLac 


5.33191 


-18.10 


13.94 


0.378 


0.194 


4.758 


3.431 


0.006 


0.006 


0.155 


0.230 


(5Cep 


5.36630 


-16.40 


15.76 


0.443 


0.293 


4.900 


3.381 


0.003 


0.003 


0.110 


0.163 


CV Mon 


5.37879 


19.28 


15.15 


0.437 


0.212 


4.672 


3.219 


0.005 


0.005 


0.152 


0.225 


XLac 


5.44453 


-25.20 


11.18 


0.161 


0.022 


4.216 


2.001 


0.013 


0.013 


0.233 


0.344 


VCen 


5.49392 


-23.60 


14.78 


0.483 


0.190 


4.531 


2.691 


0.035 


0.035 


0.418 


0.583 


YSgi- 


5.77335 


-1.57 


14.68 


0.442 


0.237 


4.851 


3.593 


0.004 


0.004 


0.129 


0.192 


CS Vel 


5.90474 


26.85 


14.93 


0.388 


0.171 


4.802 


3.522 


0.013 


0.013 


0.237 


0.346 


FMAql 


6.11423 


-5.01 


17.22 


0.403 


0.198 


4.822 


3.821 


0.006 


0.006 


0.159 


0.234 


XX Sgr 


6.42414 


12.65 


15.44 


0.492 


0.205 


4.949 


3.904 


0.002 


0.002 


0.092 


0.136 


BE Sgr 


6.63714 


8.40 


13.50 


0.466 


0.152 


5.047 


4.346 


0.001 


0.001 


0.076 


0.113 


VCar 


6.69668 


12.60 


11.75 


0.508 


0.165 


5.085 


4.615 


0.071 


0.071 


0.595 


0.819 


USgr 


6.74523 


2.70 


14.67 


0.503 


0.190 


4.939 


4.005 


0.005 


0.005 


0.139 


0.206 


V496 Aql 


6.80703 


5.10 


9.15 


0.322 


0.063 


5.082 


4.625 


0.001 


0.001 


0.076 


0.114 


XSgr 


7.01283 


-12.70 


11.45 


0.395 


0.148 


4.748 


3.459 


0.001 


0.001 


0.065 


0.098 


UAql 


7.02410 


-0.11 


14.90 


0.521 


0. 1 83 


5.005 


4.279 


0.007 


0.007 


0.176 


0.261 


17 Aql 


7.17678 


-14.80 


15.89 


0.539 


0.246 


5.169 


4.062 


0.004 


0.004 


0.131 


0.194 


WSgr 


7.59503 


-27.70 


14.78 


0.559 


0.189 


5.136 


3.921 


0.005 


0.005 


0.146 


0.217 


UVul 


7.99069 


0.24 


12.43 


0.478 


0.104 


5.321 


3.970 


0.013 


0.013 


0.236 


0.349 


SSge 


8.38209 


-13.40 


14.51 


0.544 


0.115 


5.358 


4.086 


0.002 


0.002 


0.098 


0.146 


GHLup 


9.27787 


-16.10 


4.38 


0.146 


0.000 


5.823 


0.000 


0.183 


0.183 


0.870 


-1.634 


FN Aql 


9.48151 


12.67 


14.60 


0.163 


0.156 


0.566 


5.571 


0.005 


0.005 


0.137 


0.232 


YZSgr 


9.55369 


22.88 


13.74 


0.444 


0.122 


0.016 


4.347 


0.002 


0.002 


0.083 


0.127 


SMus 


9.66007 


-0.14 


12.85 


0.479 


0.137 


5.820 


-2.260 


0.011 


0.011 


0.209 


0.338 


SNor 


9.75424 


5.54 


14.31 


0.332 


0.161 


0.255 


-1.444 


0.007 


0.007 


0.168 


0.373 


j3Dor 


9.84308 


8.87 


14.01 


0.313 


0.192 


0.337 


4.914 


0.017 


0.017 


0.263 


0.626 


f Gem 


10.15073 


7.44 


11.81 


0.243 


0.137 


0.459 


5.046 


0.001 


0.001 


0.066 


0.099 


TWNor 


10.78618 


-56.60 


18.78 


0.193 


0.218 


0.685 


5.181 


0.013 


0.013 


0.233 


1.316 


ZLac 


10.88564 


0.09 


18.43 


0.180 


0.172 


0.209 


5.059 


0.008 


0.008 


0.181 


0.269 


XXCen 


10.95337 


-16.20 


14.23 


0.340 


0.140 


0.308 


5.238 


0.051 


0.051 


0.474 


0.689 


V340 Nor 


11.28700 


-40.30 


8.58 


0.265 


0.053 


0.206 


4.441 


0.016 


0.016 


0.413 


0.391 




11.63641 


-16.90 


20.12 


0.200 


0.205 


0.049 


5.724 


0.044 


0.044 


0.462 


69.984 


UNor 


12.64371 


-22.20 


16.88 


0.135 


0.087 


0.206 


4.921 


0.044 


0.044 


0.448 


0.646 


SU Cru 


12.84760 


-30.30 


11.06 


0.009 


0.096 


3.232 


5.500 


0.067 


0.067 


1.272 


0.792 


BNPup 


13.67310 


67.57 


24.86 


0.206 


0.021 


4.495 


1.650 


0.001 


0.001 


0.069 


0.103 


TTAql 


13.75496 


3.83 


21.93 


0.139 


0.140 


4.517 


6.262 


0.001 


0.001 


0.049 


0.364 


LSPup 


14.14640 


81.50 


22.64 


0.069 


0.126 


6.180 


5.978 


0.001 


0.001 


0.075 


0.112 


VWCen 


15.03618 


-30.60 


22.72 


0.213 


0.049 


4.674 


0.357 


0.038 


0.038 


0.579 


0.603 


XCyg 


16.38633 


7.85 


22.18 


0.228 


0.058 


4.559 


0.890 


0.004 


0.004 


0.135 


0.209 


Y Oph 


17.12633 


-7.12 


7.75 


0.114 


0.034 


5.503 


4.257 


0.024 


0.024 


0.313 


0.466 


SZ Aql 


17.14071 


6.45 


25.38 


0.263 


0.039 


4.563 


1.735 


0.003 


0.003 


0.112 


0.167 


CT Car 


18.05765 


111.10 


23.42 


0.307 


0.139 


4.576 


2.444 


0.015 


0.015 


0.275 


0.389 


VY Car 


18.90728 


1.86 


22.22 


0.272 


0.112 


4.518 


2.374 


0.029 


0.029 


0.417 


0.514 


RUSct 


19.70062 


-5.16 


21.14 


0.244 


0.014 


4.753 


3.027 


0.009 


0.009 


0.236 


0.330 


RYSco 


20.32014 


-18.00 


15.85 


0.050 


0.070 


4.077 


4.236 


0.045 


0.045 


0.643 


0.672 


RZVel 


20.39690 


24.81 


22.96 


0.206 


0.172 


4.564 


1.875 


0.023 


0.023 


0.446 


0.940 


WZSgr 


21.84960 


-17.40 


23.05 


0.323 


0.138 


4.542 


2.606 


0.003 


0.003 


0.109 


0.160 


WZCar 


23.01320 


-14.40 


22.78 


0.341 


0.147 


4.468 


2.546 


0.033 


0.033 


0.598 


0.779 


VZPup 


23.17100 


63.52 


23.00 


0.303 


0.121 


4.446 


2.450 


0.001 


0.001 


0.068 


0.088 


SW Vel 


23.44313 


23.75 


26.88 


0.378 


0.200 


4.645 


2.962 


0.009 


0.009 


0.222 


0.304 


XPup 


25.96100 


72.02 


24.12 


0.363 


0.143 


4.565 


2.891 


0.001 


0.001 


0.081 


3.168 


TMon 


27.03569 


20.67 


22.57 


0.345 


0.202 


4.580 


2.983 


0.005 


0.005 


0.152 


0.220 


RY Vel 


28.13117 


-10.60 


16.18 


0.267 


0.127 


4.333 


2.464 


0.039 


0.039 


0.600 


0.602 


KQSco 


28.69580 


-30.30 


13.84 


0.940 


0.611 


3.943 


4.083 


0.038 


0.038 


0.628 


0.616 


AQPup 


30.10400 


61.02 


23.91 


0.276 


0.062 


4.733 


3.298 


0.001 


0.001 


0.066 


0.099 


KNCen 


34.02964 


-37.20 


21.71 


0.295 


0.341 


3.861 


2.323 


0.036 


0.036 


0.569 


0.655 


^ Car 


35.54804 


2.28 


16.77 


0.296 


0.150 


4.639 


2.949 


0.005 


0.005 


0.153 


0.225 


UCar 


38.81233 


0.21 


21.00 


0.347 


0.116 


4.494 


2.531 


0.022 


0.022 


0.374 


0.503 


RS Pup 


41.47446 


25.91 


19.81 


0.444 


0.192 


4.615 


3.103 


0.008 


0.008 


0.185 


0.267 


SV Vul 


44.98561 


-1.57 


19.54 


0.408 


0.197 


4.495 


3.108 


0.003 


0.003 


0.115 


0.167 


GYSge 


51.56031 


15.96 


13.15 


0.301 


0.104 


4.766 


3.118 


0.014 


0.014 


0.250 


0.362 


S Vul 


68.77426 


1.52 


13.26 


0.335 


0.208 


4.742 


3.323 


0.008 


0.008 


0.181 


0.268 
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Table 6. The derived quantities for different adopted /^-factor relations, p - /3p + Up log P, for the LMC (LMC), and Milky Way 
(MW) samples. Each column corresponds to a different p-factor relation where the slope and zero point are given in the first two 
rows. The PL relations are of the form M„ = a™ x (log(P) - 1.0) + where the index m refers to the photometric band. 



a,, 
Pp 


-0.03 
1.39 


-0.08 
1.31 


-0.08 
1.455 


-0.186 
1.550 


Parameter 




(m - M)o(LMC) 


18.50 ±0.04 


18.26 ± 0.04 


18.50 ± 0.04 


18.45 ± 0.04 


(m - M)o(LMC) slope 
M (pc) 

Mjdavg 


0.31 ±0.10 
-5 ± 11 
-0.02 ± 0.02 


0.22 ±0.10 

-40 ± 9 
-0.11 ±0.02 


0.24 ±0.10 

± 10 
0.00 ± 0.02 


0.00 ±0.10 

± 7 
0.00 ± 0.02 


bK 


-3.58 ± 0.09 
-5.65 ± 0.03 


-3.49 ± 0.08 
-5.45 ± 0.03 


-3.50 ± 0.09 
-5.67 ± 0.03 


-3.33 ± 0.09 
-5.66 ± 0.03 


civ 
bv 


-2.92 ±0.10 
-3.95 ± 0.03 


-2.83 ±0.10 
-3.73 ± 0.03 


-2.84 ±0.10 
-3.97 ± 0.03 


-2.67 ±0.10 
-3.96 ± 0.03 


ClWvi 

bwvi 


-3.51 ±0.11 
-5.95 ± 0.04 


-3.42 ±0.12 
-5.74 ± 0.04 


-3.43 ±0.11 
-5.97 ± 0.04 


-3.26 ±0.12 
-5.96 ± 0.04 


bwjk 


-3.69 ±0.10 
-5.95 ± 0.03 


-3.60 ±0.10 
-5.74 ± 0.03 


-3.61 ± 0.09 
-5.98 ± 0.03 


-3.44 ± 0.09 
-5.96 ± 0.03 
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Table 7. The distances and intensity-averaged absolute magnitudes for the complete sample of Milky Way Cepheids based on 
the IRSB method as caUbrated in this paper. The cr values are the nominal values returned by the bi-sector fitting algorithm. The 
Wesenheit indices based on the (V, /) and {J, K) photometry are given as well as the adopted reddening. The last column gives the 
adopted phase shift, A0, between spectroscopic and photometric data. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


ID 


logP 


d 


a-(d) 


(m - M)o 




Mb 


Mv 


M, 


M] 


M„ 


Mk 


Wvi 


WjK 


E{B - V) 


hip 






(pc) 


(pc) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 




SU Cas 


0.289884 


418 


12 


8.10 


0.06 


-2.53 


-2.97 


-3.52 


-3.82 


-4.05 


-4.08 


-4.37 


-4.25 


0.259 


0.000 


DTCyg 


0.397804 


621 


30 


8.97 


0.11 


-2.82 


-3.32 


-3.87 


-4.25 


-4.51 


-4.54 


-4.71 


-4.75 


0.040 


0.030 


EVSct 


0.490098 


1766 


82 


11.24 


0.10 


-2.72 


-3.21 


-3.90 


-4.19 


-4.45 


-4.46 


-4.97 


-AM 


0.655 


0.045 


SZTau 


0.498166 


558 


6 


8.73 


0.02 


-2.62 


-3.16 


-3.80 


-4.18 


-4.49 


-4.53 


-4.78 


-4.77 


0.295 


-0.005 


SS Set 


0.564814 


1000 


33 


10.00 


0.07 


-2.21 


-2.83 


-3.55 


-3.92 


-4.24 


-4.27 


-4.65 


-4.51 


0.325 


-0.035 


RT Aur 


0.571489 


389 


4 


7.95 


0.02 


-2.16 


-2.69 


-3.30 


-3.75 


-4.01 


-4.05 


-4.23 


-4.26 


0.059 


0.005 


QZNor 


0.578244 


1346 


38 


10.64 


0.06 


-2.00 


-2.62 




-3.76 


-4.07 


-4.12 




-4.38 


0.263 


0.020 


SU Cyg 


0.584952 


909 


21 


9.79 


0.05 


-2.77 


-3.24 


-3.79 


-4.17 


-4.45 


-4.50 


-4.64 


-4.73 


0.098 


0.000 


BFOph 


0.609329 


708 


10 


9.25 


0.03 


-2.04 


-2.67 


-3.35 


-3.77 


-4.10 


-4.16 


-4.39 


-4.43 


0.235 


0.030 


YLac 


0.635863 


2430 


35 


11.93 


0.03 


-2.97 


-3.48 


-4.08 


-4.43 


-4.74 


-4.79 


-4.99 


-5.04 


0.217 


-0.005 


TVul 


0.646934 


543 


5 


8.67 


0.02 


-2.54 


-3.12 


-3.74 


-4.13 


-4.44 


-4.50 


-4.70 


-4.75 


0.060 


0.020 


FFAql 


0.650390 


369 


11 


7.84 


0.06 


-2.54 


-3.10 


-3.72 


-4.10 


-4.38 


-4.42 


-4.68 


-4.63 


0.196 


-0.030 


T Vel 


0.666501 


1002 


11 


10.00 


0.02 


-2.28 


-2.91 




-4.05 


-4.40 


-4.47 




-4.76 


0.289 


0.000 


VZCyg 


0.687034 


1897 


35 


11.39 


0.04 


-2.68 


-3.29 


-3.98 


-4.36 


-4.68 


-4.74 


-5.04 


-5.01 


0.266 


0.035 


V350 Sgr 


0.712165 


990 


23 


9.98 


0.05 


-2.86 


-3.47 


-4.15 


-4.55 


-4.89 


-4.93 


-5.20 


-5.19 


0.299 


0.000 


BGLac 


0.726883 


1776 


26 


11.25 


0.03 


-2.68 


-3.33 


-4.05 


-4.42 


-4.76 


-4.82 


-5.16 


-5.09 


0.300 


0.005 


5Cep 


0.729678 


267 


5 


7.13 


0.04 


-2.83 


-3.41 


-4.07 


-4.45 


-4.78 


-4.84 


-5.08 


-5.11 


0.072 


0.000 


CVMon 


0.730685 


1547 


24 


10.95 


0.03 


-2.40 


-2.97 


-3.75 


-4.22 


-4.58 


-4.65 


-4.94 


-4.94 


0.722 


0.020 


XLac 


0.735997 


1852 


38 


11.34 


0.04 


-3.45 


-4.02 


-4.68 


-5.05 


-5.38 


-5.42 


-5.70 


-5.68 


0.336 


0.040 


VCen 


0.739882 


643 


20 


9.04 


0.07 


-2.58 


-3.16 


-3.82 


-4.25 


-4.58 


-4.65 


-4.85 


-4.93 


0.294 


0.000 


YSgr 


0.761428 


437 


15 


8.20 


0.07 


-2.40 


-3.07 


-3.78 


-4.25 


-4.61 


-4.67 


-4.87 


-4.95 


0.188 


-0.030 


CS Vel 


0.771201 


3184 


108 


12.52 


0.07 


-2.59 


-3.19 




-4.37 


-4.71 


-4.78 




-5.07 


0.737 


0.015 


FM Aql 


0.786342 


1182 


28 


10.36 


0.05 


-3.31 


-4.00 


-4.80 


-5.12 


-5.47 


-5.52 


-6.05 


-5.79 


0.589 


-0.010 


XX Sgr 


0.807815 


1042 


42 


10.09 


0.09 


-2.33 


-2.92 


-3.65 


-4.11 


-4.44 


-4.47 


-4.77 


-4.73 


0.521 


-0.025 


BB Sgr 


0.821971 


868 


13 


9.69 


0.03 


-2.96 


-3.66 


-4.42 


-4.86 


-5.22 


-5.29 


-5.60 


-5.59 


0.281 


0.000 


VCar 


0.825860 


869 


16 


9.69 


0.04 


-2.16 


-2.86 




-4.04 


-4.40 


-4.47 




-4.77 


0.166 


0.000 


USgr 


0.828997 


579 


6 


8.81 


0.02 


-2.73 


-3.42 


-4.18 


-4.61 


-4.96 


-5.02 


-5.34 


-5.30 


0.402 


0.000 


V496 Aql 


0.832958 


829 


26 


9.59 


0.07 


-2.37 


-3.12 




-4.32 


-4.68 


-4.70 




-4.96 


0.397 


-0.055 


XSgr 


0.845893 


322 


5 


7.54 


0.04 


-3.23 


-3.74 


-4.35 


-4.75 


-5.06 


-5.11 


-5.28 


-5.36 


0.237 


0.000 


UAql 


0.846591 


592 


19 


8.86 


0.07 


-2.91 


-3.58 


-4.31 


-4.72 


-5.06 


-5.11 


-5.44 


-5.38 


0.360 


0.015 


J/Aql 


0.855930 


255 


4 


7.03 


0.09 


-2.89 


-3.55 


-4.27 


-4.70 


-5.04 


-5.10 


-5.38 


-5.38 


0.129 


0.000 


W Sgr 


0.880529 


217 


16 


6.68 


0.16 


-1.73 


-2.37 


-3.03 


-3.47 


-3.80 


-3.88 


-4.05 


-4.17 


0.108 


0.040 


U Vul 


0.902584 


658 


12 


9.09 


0.04 


-3.22 


-3.92 


-4.73 


-5.03 


-5.34 


-5.37 


-6.00 


-5.60 


0.603 


0.005 


S Sge 


0.923352 


671 


16 


9.13 


0.05 


-3.13 


-3.84 


-4.55 


-5.00 


-5.34 


-5.41 


-5.66 


-5.70 


0.099 


0.000 


GHLup 


0.967448 


1307 


23 


10.58 


0.04 


-3.21 


-4.07 


-4.90 


-5.42 


-5.82 


-5.90 


-6.18 


-6.23 


0.347 


-0.040 


FN Aql 


0.976878 


1175 


20 


10.35 


0.04 


-2.81 


-3.53 


-4.36 


-4.75 


-5.12 


-5.17 


-5.63 


-5.46 


0.486 


-0.030 


YZSgr 


0.980171 


1137 


13 


10.28 


0.03 


-3.10 


-3.84 




-5.07 


-5.44 


-5.49 




-5.77 


0.281 


-0.010 


SMus 


0.984980 


858 


17 


9.67 


0.04 


-3.30 


-3.98 




-5.25 


-5.62 


-5.71 




-6.02 


0.140 


0.010 


SNor 


0.989194 


950 


11 


9.89 


0.02 


-3.27 


-4.04 


-4.81 


-5.33 


-5.72 


-5.79 


-6.01 


-6.12 


0.179 


0.000 


ySDor 


0.993131 


327 


5 


7.57 


0.03 


-3.25 


-4.00 


-4.73 


-5.18 


-5.57 


-5.63 


-5.85 


-5.94 


0.051 


0.000 


^ Gem 


1.006497 


386 


9 


7.93 


0.05 


-3.29 


-4.13 




-5.33 


-5.74 


-5.81 




-6.13 


0.014 


0.050 
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Fig. 11. The difference in derived distance modulus A(m - M) - 
(m - M)zAMS - (ni - M)irsb- The open symbol shows the star 
SU Cas which has been disregarded in the comparison. 



4. Discussion 

The slopes of the PL relations which we find depend directly on 
the adopted /^-factor relation and as we have shown the revised, 
empirically-determined relation is necessary to give pulsation 
period-independent distances to LMC Cepheids. In Paper II we 
further show that the adopted p-factor relation also reproduces 
the slopes of the apparent magnitude versus log(P) relations both 
in the near-IR and in the optical bands. The revised /^-factor re- 
lation confirms our earlier findings (Gieren et al. 2005 ) based on 
a much smaller sample of LMC Cepheids. Still, it is at odds with 
the recent careful theoretical study by Nardetto et al. (120091 ) . We 
do not have a ready explanation for this disagreement but sus- 
pect that it is either due to some missing physics in the theo- 
retical approach or some previously undetected period depen- 
dence in the IRSB method, which we have eliminated through 
the parametrization chosen here. If the problem does not orig- 
inate from the pulsational velocities, i.e. with the p-factor, it 
would have to originate with the surface-brightness calibration. 
The calibration is based on mostly short period Cepheids with 
only a single long period (35 days) Cepheid. However, Kervella 
et al. (I2004bl l show a quite convincing comparison of surface- 
brightness relations for short and long period stars and the fact 
that the relation agrees well with the relation for static stars is 
also suggestive that the cause for the steep slope is not buried 
here. The remaining problem seems to be with the short period 
stars where the revised p-factor relation leads to values which get 
close to L5, suggesting a uniform disk and no limb-darkening. 
However, the short period stars are the ones which are most nu- 
merous in the sample used by Kervella et al. (I2004bl) and thus the 
ones which have the best empirical constraints. Still we intend 
to make direct comparisons with interferometric studies as done 
by Kervella et al. (2004c ) for { Car in an attempt to understand 
the reason for the effect which we see. For our main purpose of 
investigating the effect of metallicity on the Cepheid PL relation 
by comparing Milky Way and Magellanic Cloud Cepheids we 
are working in a purely differential way so the actual p-factor 
relation cancels out as long as it is the same for both samples of 
stars, i.e. metal independent, which from theoretical considera- 
tions seems to be a good assumption (Nardetto et al. l201 II ). 



Groenewegen (120071 ) argues for a constant p-factor relation 
based on a comparison with a subset of the HST parallax stars 
but the scatter in his Fig. 2 is very large. With our new high pre- 
cision radial velocity data for three of the HST parallax stars, 
we confirm our p-factor relation (see Fig|7]) through excellent 
agreement with the full set of Benedict et al. (120071 ) results on a 
star by star basis. We also base our p-factor relation on a much 
larger sample of LMC Cepheids. This sample has a large number 
of stars with pulsation periods significantly longer than ten days, 
thus forming a much firmer basis for constraining the p-factor 
relation. In a following paper, Groenewegen (120081 ) discusses 
the use of the Nardetto et al. (120071 ) relation and finds no signif- 
icant difference to a constant value. At the same time he finds 
period-luminosity relations for Milky Way Cepheids which are 
very similar to the relations presented in the previous section. 
He finds in the /T-band a relation with a slope of -3.38 + 0.08, 
in good agreement with our value of -3.33 + 0.09, similarly he 
finds a slope in the V-band of -2.60 + 0.09 whereas we find a 
value of -2.67 + 0.10. It thus seems as if, in spite of the fact 
that we apply the same method, there are significant differences 
in the implementation of the technique which might affect the 
results. 

Recently Molinaro et al. (1201 11), based on the CORS variant 
of the Baade-Wesselink method and using Walraven photometry 
for 26 galactic Cepheids and a constant p-factor of L27, found a 
PL relation of = -2.78(+0.11) log(P) - 1.42(+0.11). Again 
the slope is in good agreement with our relation whereas the zero 
point at a period of ten days differs by 0.23 mag, our value being 
fainter. 

There are presently very few alternative routes to delineat- 
ing the Milky Way PL relation apart from the Baade-Wesselink 
type methods. There is of course the recent direct parallaxes to 
ten Cepheids with the HST fine guidance sensors by Benedict et 
al. (2007 ), but this constitutes a modest sample of stars for a PL 
relation. The classical approach is the zero-age main sequence 
(ZAMS) fitting to OB associations and open clusters contain- 
ing Cepheids (see e.g. Feast and Walker (119871 ) and references 
therein). Turner (2010) rederived the Milky Way PL relation 
based on the ZAMS fitting to OB associations and open clus- 
ters containing Cepheids finding My = -2.78(+0.12) log(P) - 
1.29(+0.10). In Figdllwe compare the ZAMS fitting based dis- 
tance moduli with the IRSB based moduli for the stars in com- 
mon, and we find very good agreement with no significant period 
dependence. The unweighted mean difference is 0.12 + 0.06 mag 
so we do find a slight zero point offset. We have excluded the star 
SU Cas in the comparison as it is an outlier, even if the IRSB fit 
appears well defined and does not indicate any obvious problem 
with the data. 

Benedict et al. (20071) found a slope of -2.43 + 0.12 in the V 
band from the HST parallax measurements, a value which is only 
slightly shallower than our value of -2.67 + 0.10 and certainly 
not steeper than our value. 

We argue that most recent investigations agree to within the 
errors with the slope of our relation and they disagree with the 
earlier findings of Sandage et al. (120041 ) and Storm et al. (120041 ). 
that the Milky Way PL relations are significantly steeper than the 
LMC relations. In fact, in Paper II we find that if anything the 
optical Milky Way PL relations might be slightly shallower than 
the LMC relations. 

In Paper II we find that both the slopes and the zero-points 
of the near-IR PL relations are insensitive to metallicity. In ad- 
dition the /T-band PL relation is very insensitive to reddening 
making this relation our preferred standard candle. We argue 
that the best calibration of this relation is the combined /iT-band 
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PL relation given in Tab|9] with the small metallicity effect of 
J - -0.10 + 0.10 mag dex"' also given in that table. We note 
that for most extra-galactic Cepheid samples the metallicity is 
close to the range from LMC to solar and the failure to correct 
for the metallicity effect leads to systematic errors of the order 
of only 0.02 mag in the distance modulus. 

In the optical bands, the LMC and Milky Way slopes are 
less in agreement differing by up to 0.2 mag dex"' as shown in 
Paper II. However in the compilation by Bono et al. (1201 Oi l the 
slopes in the V and / bands for extra-galactic samples show a 
large spread of the order 1 mag dex"' , much larger than our ob- 
served difference between the LMC and MiUcy Way samples. 
The slope variations in that paper do not seem to be strongly 
correlated with metallicity variations so from that point of view 
our combined SMC, LMC, MW relation also provides the better 
reference relation as it is based on more stars. 



5. Conclusions and Summary 

We have obtained new, accurate, radial velocity curves for four- 
teen Milky Way Cepheids including three Cepheids with direct 
parallax measures from Benedict et al. (2007), expanding the 
sample of Milky Way fundamental mode Cepheids to which we 
can determine precise IRSB distances to a total of 70 stars. 

We have empirically redetermined the p-factor relation, 
which converts the observed radial velocities into pulsation ve- 
locities needed for the IRSB method, using two fundamental 
physical constraints. The first constraint is that the distance to 
LMC Cepheids should be independent of their pulsation peri- 
ods, and the second constraint is that on average we should 
reproduce the distances to the Cepheids with parallaxes from 
Benedict et al. C2007 ). We find quite a steep relation, p - 
1.550(+0.04) - 0.186(+0.06)log(P) which is not easily recon- 
ciled with recent theoretical work (e.g Nardetto et al. 120091 ). 
However, this revised relation gives rise to PL relations which 
are in excellent agreement with other independent determina- 
tions both for the Milky Way, as shown in the present paper, and 
for the LMC, as shown in Paper II. 

Using the revised p-factor relation we have determined pre- 
cise PL relations in the V, /, /, & K bands, as well as the 
Wesenheit indices Wy/ , & Wjk for these Milky Way stars. These 
relations can be used for distance determination to other galaxies 
with solar abundance. 

In Paper II we compared these relations to similar rela- 
tions for a sample of LMC Cepheids and we found that the 
effect of metallicity on the slopes is negligible in the near-IR 
and small, possibly consistent with zero, in the optical bands as 
well. Including also a sample of SMC Cepheids we find that the 
zero points of the PL relations depend on metallicity to a vary- 
ing degree, but in most bands the effect is small, of the order 
-0.10 + 0.10 mag dex"' , which is consistent with a zero effect. 
Consequently we argue that it is warranted to combine the three 
samples of Cepheids giving us a total sample of 1 1 1 Cepheids 
with IRSB distances which can be used to delineate accurate ab- 
solute, universal PL relations. 

Our best standard candle is the /T-band PL relation as it is not 
only insensitive to reddening and shows a low intrinsic disper- 
sion, but it also exhibits no metallicity dependence on the slope 
and only a weak dependence on the zero point consistent with 
a null effect. The TiT-band relation based on the full sample of 
Cepheids presented here is Mk = -3.30(+0.06)[log(P) - 1.0] - 
5.65(+0.02). The combined relation in the optical Wv/ index is 
Wyi = -3.32(+0.08)[log(P) - 1.0] - 5.92(+0.03). 
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